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What is commonly called the reheat temperature, Trh, is not the max- 
imum temperature obtained after inflation. The maximum temperature is, 
in fact, much larger than Trh- As an application of this we consider the pro- 
duction of massive stable dark-matter particles of mass Mx during reheating, 
and show that their abundance is suppressed as a power of T RH /M X rather 
than exp(— Mx/Trh). We find that particles of mass as large as 2 x 10 3 times 
the reheat temperature may be produced in interesting abundance. In ad- 
dition to dark matter, our analysis is relevant for baryogenesis if the baryon 
asymmetry is produced by the baryon (or lepton) number violating decays of 
superheavy bosons, and also for relic ultra-high energy cosmic rays if decays 
of superheavy particles are responsible for the highest energy cosmic rays. 
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I. INTRODUCTION 



At the end of inflation [1] the energy density of the universe is locked up in a combi- 
nation of kinetic energy and potential energy of the inflaton field, with the bulk of the 
inflaton energy density in the zero-momentum mode of the field. Thus, the universe at the 
end of inflation is in a cold, low-entropy state with few degrees of freedom, very much un- 
like the present hot, high-entropy universe. After inflation the frozen inflaton-dominated 
universe must somehow be defrosted and become a high-entropy radiation-dominated 
universe. 

The process by which the inflaton energy density is converted to radiation is known 
as "reheating" ||. The possible role of nonlinear dynamics leading to explosive particle 
production has recently received a lot of attention. This process, known as "preheat- 
ing" || may convert a fair fraction of the inflaton energy density into other degrees of 
freedom, with extremely interesting cosmological effects such as symmetry restoration, 
baryogenesis, or production of dark matter. But the efficiency of preheating is very sen- 
sitive to the model and the model parameters. In some models the process is inefficient; 
in some models it is not operative at all. Even if preheating is relatively efficient, it is 
unlikely to remove all of the energy density of in the inflaton field. It is likely that the 
slow decay of the inflaton field is necessary to extract the remaining inflaton field energy. 

The simplest way to envision this process is if the comoving energy density in the zero 
mode of the inflaton decays into normal particles, which then scatter and thermalize to 
form a thermal background. It is usually assumed that the decay width of this process 
is the same as the decay width of a free inflaton field. 

There are two reasons to suspect that the inflaton decay width might be small. The 
requisite flatness of the inflaton potential suggests a weak coupling of the inflaton field 
to other fields since the potential is renormalized by the inflaton coupling to other fields 
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|T||. However, this restriction may be evaded in supersymmetric theories where the non 



renormalization theorem ensures a cancelation between fields and their superpartners. 
A second reason to suspect weak coupling is that in local supersymmetric theories grav- 
itinos are produced during reheating. Unless reheating is delayed, gravitinos will be 
overproduced, leading to a large undesired entropy production when they decay after 
big-bang nucleosynthesis @] . 

Of particular interest is a quantity known as the reheat temperature, denoted as 
Trh- The reheat temperature is calculated by assuming an instantaneous conversion 
of the energy density in the inflaton field into radiation when the decay width of the 
inflaton energy, F^, is equal to H, the expansion rate of the universe. 

The reheat temperature is calculated quite easily 0. After inflation the inflaton field 
executes coherent oscillations about the minimum of the potential. Averaged over several 
oscillations, the coherent oscillation energy density redshifts as matter: p^ oc a~ 3 , where 
a is the Robertson-Walker scale factor. If we denote as pi and aj the total inflaton 
energy density and the scale factor at the initiation of coherent oscillations, then the 
Hubble expansion rate as a function of a is (Mpi is the Planck mass) 



Equating H(a) and leads to an expression for aj/a. Now if we assume that all 
available coherent energy density is instantaneously converted into radiation at this value 
of ai/a, we can define the reheat temperature by setting the coherent energy density, 
P(f> — Pi(ai/ a ) 3 , equal to the radiation energy density, p^ = (n 2 / 30) g*T^ H , where g* is 
the effective number of relativistic degrees of freedom at temperature T RH . The result is 




(1) 





The limit from gravitino overproduction is Trh < 10 9 to 10 10 GeV. 
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The reheat temperature is best regarded as the temperature below which the universe 
expands as a radiation-dominated universe, with the scale factor decreasing as g* 1 T~ x . 
In this regard it has a limited meaning 0, |5|. For instance, T R h should not be used as 
the maximum temperature obtained by the universe during reheating. The maximum 
temperature is, in fact, much larger than Trh- One implication of this is that it is 
incorrect to assume that the maximum abundance of a massive particle species produced 
after inflation is suppressed by a factor of exp(— M/T R h). 

In this paper we illustrate this effect by calculating the abundance of a massive parti- 
cle species produced in reheating. We show that particles of mass much greater than the 
eventual "reheat" temperature T RH may be created by the thermalized decay products 
of the inflaton. As an example, we demonstrate that a stable particle species X of mass 
Mx would be produced in the reheating process in sufficient abundance that its contri- 
bution to closure density today is approximately M x (a\v\)(g*/200)~ 3 / 2 (2000T R H / M x ) 7 , 
where g* is the number of effective degrees of freedom of the radiation energy density 
and (cx|v|) is the thermal average of the X annihilation cross section times the M0ller 
flux factor. Thus, particles of mass as large as 2000 times the reheat temperature may 
be produced in interesting abundance. 

Other applications of the effect include production of massive Higgs bosons which 
could decay and produce the baryon asymmetry, or massive particles that could decay 
and produce high-energy cosmic rays. 

In the next section we develop a system of Boltzmann equations describing the evolu- 
tion of the energy densities of the inflaton field, radiation, and a massive particle species. 
In Section III we find analytic approximations to the system and estimate the contribu- 
tion to the present critical density from a stable, massive particle produced in reheating. 
Section IV contains some numerical results which illustrate several generic features. We 
conclude in Section V by discussing some applications of our results. The assumption 
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of local thermodynamic equilibrium for the light degrees of freedom is addressed in an 
appendix. 



II. BOLTZMANN EQUATIONS DESCRIBING REHEATING 

Let us consider a model universe with three components: inflaton field energy, p^, 
radiation energy density pr, and the energy density of a nonrelativistic particle species, 
px- We will assume that the decay rate of the inflaton field energy density is T^, with a 
branching fraction into XX of Bx, and a branching fraction 1 — Bx into light degrees of 
freedom, generically referred to as radiation. We will denote the decay width of the X 
as Fx- We will also assume that the light degrees of freedom are in local thermodynamic 
equilibrium. This is by no means guaranteed, and we will return to the question in the 
appendix. 

With the above assumptions, the Boltzmann equations describing the redshift and 
interchange in the energy density among the different components is 



p<p + ZHp^ + F^p^ = 

p R + AH p R - (1 - B x )F m - 

m x 

(a\v\) 



Px-(p E x Q ) 2 }-r x (p x -Px Q )=o 



px + SHpx - BxF^p^ + 



Px-ipT) +rx{px-P E x Q )=o, (3) 



m x 

where dot denotes time derivative As already mentioned, (cr\v\) is the thermal average 
of the X annihilation cross section times the M0ller flux factor. The equilibrium energy 
density for the X particles, p x Q , is determined by the radiation temperature, T. 

It is useful to introduce two dimensionless constants, and ax, defined in terms of 

and (a\v\) as 

= a M (a\v\) = a x M x 2 . (4) 
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For a reheat temperature much smaller than Ms, must be small. From Eq. (0), the 



reheat temperature in terms of ax and Mx is Trh — a l J 2 J M^Mpi. For = 10 13 GeV, 
a^ must be smaller than of order 10~ 13 . On the other hand, ax may be as large as of 
order unity, or it may be small also. 

In what follows we will make the simplifying assumption that Bx = 0. Since we are 
interested in a stable particle relic, we will assume that Tx = 0. It is also convenient 
to work with dimensionless quantities that can absorb the effect of expansion of the 
universe. This may be accomplished with the definitions 

$ = P4> M^a 3 ; R = p R a 4 ; X = pxM^a 3 . (5) 

It is also convenient to use the scale factor, rather than time, for the independent variable, 
so we define a variable x = aM^. With this choice the system of equations can be written 
as (prime denotes d/dx) 

$' = -a , x $ 

' 2 =$ +C 2 ; = (X 2 ~X 2 

x~ 2 



R' = a - ' $ + c 2 ; (x 2 - x 2 EO 

X' = -c 3 /" 2 (X 2 -X 2 ) . (6) 



The constants ci, c 2 , and c 3 are given by 

~-Mn M^ax M, : , frj . 

V o7r M x M x 

Xeq is the equilibrium value of X, given in terms of the temperature T as (assuming a 

single degree of freedom for the X species) 

M 3 / 1 \ 3 / 2 / T \ 3/2 
Xe Q = ^{-) x 3 —-) eM-Mx/T). (8) 



M 3 \2nJ \M X . 

The temperature depends upon R and g*, the effective number of degrees of freedom in 
the radiation: 

T ( 30 \ 1/4 i?V4 

(9) 



M x \g*n 2 J M x x 



It is straightforward to solve the system of equations in Eq. (||) with initial conditions 
at x — Xi of R(xi) = X(xi) = and = $j. It is convenient to express p${x = xi) 

in terms of the expansion rate at xi, which leads to 

3 M 2 Pl H 2 



8vr Ml Ml 



X 1 



(10) 



The numerical value of xj is irrelevant. 

Before numerically solving the system of equations, it is useful to consider the early- 
time solution for R. Here, by early time, we mean H ^> r^, i.e., before a significant 
fraction of the comoving coherent energy density is converted to radiation. At early 
times $ ~ $/, and R ~ X ~ 0, so the equation for R! becomes R' = cix 3 ^ 2 ^ 1 / 2 . Thus, 
the early time solution for R is simple to obtain: 



R ~ -c x (x^ 2 - x) /2 ) $} /2 (H » T 4 



Now we may express T in terms of R to yield the early-time solution for T: 



(11) 



T 

~Mh 



12 

n 2 g* 



1/4 



c 1/4 



$1 



X 1 



X 



Xj 



-3/2 



X 



Xj 



-,1/4 



(#»r ? 



(12) 



Thus, T has a maximum value of 



Ma, 



"(*)'"-"(S)" 



0.77a 



1/4 
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27r 3 g* 



1/4 (Iff 



(13) 



which is obtained at x/xi = (8/3) 2/,s = 1.48. It is also possible to express in terms of 
Trh and obtain 



Tmax _ qjj I 9 



'HjM Pl 



1/4 



T 2 



' V57r 3 5f* J 

For an illustration, in the simplest model of chaotic inflation Hj ~ with M, 
10 13 GeV, which leads to T MAX /T RH ~ 10 3 (200/^) 1 / 8 for Tr H = 10 9 GeV. 



(14) 



We can see from Eq. (|TT| ) that for xjxj > 1, in the early-time regime T scales as 
a~ 3 / 8 , which implies that entropy is created in the early-time regime [[J. So if one is 
producing a massive particle during reheating it is necessary to take into account the 
fact that the maximum temperature is greater than Trh, and that during the early-time 
evolution, T oc a -3 / 8 . 



III. PRODUCTION OF A MASSIVE, STABLE PARTICLE SPECIES 
A. Freeze out of the comoving X energy density 

In this section we develop the equation for the Xp, the final value of X, which can 
be found from the early-time behavior. 

At early times $ ~ and R ~ 0. We will here also assume that X <C Xeq- 
Numerical results confirm the validity of this approximation and show that the massive 
particles are never in chemical equilibrium (although presumably they are in kinetic 
equilibrium). The early-time equation for the development of the X energy density is 

X' = c^\-^Xl Q . (15) 

Xeq is given in terms of Mx and the temperature, which may be found from the early- 
time solution for R. 

We can integrate Eq. (|15|) by approximating it as a Gaussian integral. First we 
rearrange Eq. fll~5|) by making appropriate redefinitions. Define the quantities y and v 
by y = X/xj and v = (x/xj) 3 ^ 16 . Now, using Eq. (|12D, we can rewrite Eq. (|15|) as 

yiy) = Qj"di/ exp[-H(i/)] (16) 
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where we have defined 

3 l/2 2 l/4\ M 3/2 M 



_ 3/4 / O-'-Z-' ' \ 1V1 PI 1V1 X 

Q = a J a x I - ' 



H{v) = A(z/- 8 -z/- 64 / 3 )~ 1/4 -|ln(z/- 8 -z/- 64 / 3 ) -231m/ 



2 5 / 4 vr 3 / 4 ^ 1/4 a, 1/4 M 



To proceed with the Gaussian integral approximation, we assume z/g~ 8 ^> ^o~ 64//3 
where i/q is the solution to H'(uq) = 0. Then, we can easily solve H'(uq) = 0, finding 



uo = J 17/2A, which is the point about which we Taylor expand H(y) to quadratic order. 
Since the integrand falls to rapidly away from u = u , and since we desire the freeze 
out value for y, the limits of the integrand can be taken to ±oo. We thereby find 
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X F ^ _ 3 B a$a x ( H]M Pl \ ( \/Yl\ 

l$*y~*^![-M^)[-^) ex P (-17/2). (18) 



Using Eq. @, we rewrite this in a more transparent form 



X F _ 4.21 x 10-6 axHfMfA hq) 
x\ ~ (^/200)3/2 M*M% ' 1 ) 

Note that this approximation should be valid as long as Uq 8 ^> isq 64 ^ 3 is satisfied. If 
the condition is not satisfied, the suppression will be exponential in Mx/T. 



B. VLxh 2 in terms of Xp 

After freeze out of the comoving energy density of the stable particle, X remains con- 
stant, so px(% > xp) = X F xJ 3 M x M^(xi/x) 3 . For delayed reheating (T^ <C Hj) freeze 
out will be well before reheating. After reheating, px(x > xrh) = Px(xrh){xrh/x) 3 - 
The comoving entropy density is constant after reheating, so the radiation energy density 
scales as pn(x > xrh) = Pr(xrh)[9*(Trh) / Q^T)] 1 / 3 (xrh / ■ Using these facts, we can 
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express the present contribution of the massive particle species to the critical density in 
terms of the ratio of the energy densities at freeze out: 

n xh 2 = p x (T RH ) f g, (today) \ 1/3 _Xo_ = p x (T RH )T RH ^ > ^ 

Q R h 2 Pr{T rh ) \ g*(T RH ) ) x RH Pr{T rh ) T r '"' 

where xq is the present value of x and To = 2.37 x 10~ 13 GeV is the present temperature.0 
Today, Q R h 2 = 4.3 x 10~ 5 , and the contribution to Qh 2 from the massive particle is 

n 12 , <: ,,,i S (T S H\X F MxMl ,„,. 
= X 10 (loW ^fi^ ' (21) 

Using the expression for Xpjx\ from the previous section, we arrive at the final result 



IV. NUMERICAL RESULTS 

An example of a numerical evaluation of the complete system in Eq. (|]) is shown 
in Fig. |1| The model parameters chosen were = 10 13 GeV, a$ = 2 x 10 -13 , Mx = 
1.15 x 10 12 GeV, ax = 10~ 2 , and = 200. The expansion rate at the beginning of 
the coherent oscillation period was chosen to be Hj = M$. These parameters result in 
T RH = 10 9 GeV and Vt x h 2 = 0.3. 

Figure [l| serves to illustrate several aspects of the problem. Just as expected, the 
comoving energy density of <fr (he., a 3 p ( / ) ) remains roughly constant until ~ H, which 
for the chosen model parameters occurs around a/aj ~ 5 x 10 8 . But of course, that does 
not mean that the temperature is zero. Notice that the temperature peaks well before 

"reheating." The maximum temperature, T M ax = 10 12 GeV, is reached at a/a I slightly 

2 In this subsection we make the heretofore criticized approximation that the infiaton energy density 
scales like pressureless matter until it dumps all of its energy into radiation at the instant of "reheating." 
In this instance, however, it is an appropriate approximation, as borne out by analytic approximations 
and the numerical calculations presented in the next section. 
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Figure 1: The evolution of energy densities and T/Mx as a function of the scale factor. 
Also shown is X/X E q. 

larger than unity (in fact at a/ai = 1.48 as expected), while the reheat temperature, 
Trh = 10 9 GeV, occurs much later, around a/ai ~ 10 8 . Note that Tmax — 10 3 Trh in 
agreement with Eq. (|i~4T). 

From the numerical results we can justify one of the assumptions in deriving the 
analytical approximations. From the figure it is clear that X <ti Xeq at the epoch of 
freeze out of the comoving X number density, which occurs around a/aj ~ 10 2 . The 
rapid rise of the ratio after freeze out is simply a reflection of the fact that X is constant 
while Xeq decreases exponentially. The relevance of the ratio is the justification of the 
neglect of Xeq term in Eq. fllpp. 

A close examination of the behavior of T shows that after the sharp initial rise of the 
temperature, the temperature decreases as a~ 3 / 8 [as follows from Eq. (JfJD] until H ~ r^, 
and thereafter T oc a -1 as expected for the radiation-dominated era. 
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For the choices of M^, a^, g*, and ax used for the model illustrated in Fig. [L], 
Qxh 2 = 0.3 for Mx = 1.15 x 10 12 GeV, in excellent agreement with the mass predicted 
by using Eq. (p^). 

V. CONCLUSIONS 

Let us now analyze the implications of our findings for the GUT baryogenesis scenario, 
where the baryon asymmetry is produced by the baryon (or lepton) number violating 
decays of superheavy bosons ||. At the end of inflation the Universe does not contain 
any matter and, even more important, it is perfectly baryon symmetric — there is no dom- 
inance of matter over antimatter. This means that GUT baryogenesis may be operative 
only if the supermassive GUT bosons are regenerated during the stage of thermalization 
of the decay products of the inflaton field[]. A naive estimate would lead to conclude that 
that the maximum number density of a massive particle species X produced after infla- 
tion is suppressed by a factor of (M x /T RH ) 3 ^ 2 exp(— Mx/Trh) with respect to the photon 
number density. For such a reason, it is commonly believed that GUT baryogenesis is 
incompatible with models of inflation where the reheating temperature is much smaller 
than the GUT scale and, in general, than the mass of the X particles, Trh ^ Mx- In 
fact, we have seen that the reheat temperature has a limited meaning and should not 
be used as the maximum temperature obtained by the universe during reheating. The 



maximum temperature, Eq. (14]), is much larger than Trh, and particles of mass much 
greater than the eventual reheating temperature Trh may be created by the thermalized 
decay products of the inflaton without any exponential suppression factor. Indeed, the 
number density nx of particles X after freeze out and reheating may be easily inferred 



3 In the case in which reheating is anticipated by a period of preheating, superheavy bosons may be 
produced by the phenomenon of parametric resonance M . 
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from Eqs. (pj) and and reads 



^3xl0-(™f 2 f^)7^). (23) 



In theories that invoke supersymmetry to preserve the flatness of the inflaton poten- 
tial, the slow decay rate of the gravitinos, the superpartners of the gravitons, is a source 
of the cosmological problems because the decay products of the gravitino will destroy 
the 4 He and D nuclei by photodissociation, and in the process destroy the successful 
nucleosynthesis predictions. The most stringent bound comes from the resulting over- 
production of D + 3 He, which would require that the gravitino abundance is smaller 
than about 10~ 10 relative to the entropy density at the time of reheating after infla- 
tion. This translates into an upper bound on the reheating temperature after inflation, 

It is easy to check that for such small values of Trh the ratio in Eq. (|23|) is always much 
larger than the equilibrium value, nf Q /n 7 = (M x /2T RH ) 3/2 (n 1/2 /£(3)) exp(-M x /T RH ). 
This result is crucial for the out-of-equilibrium decay scenarios of baryogenesis. For 
instance, in theories where B — L is a spontaneously broken local symmetry, as suggested 
by 5*0(10) unification, the cosmological baryon asymmetry can be generated by the out- 
of-equilibrium decay of the lightest heavy Majorana right-handed neutrino iVf |§, whose 
typical mass is about 10 10 GeV. For reheat temperatures of the order of 10 9 GeV, the 
number density of the right-handed neutrino is about 3 x 10 -2 n 7 and one can estimate 
the final baryon number to be of the order of B ~ {nNi/ n "t){ e / 9*) ~ 10~ 4 e, where e 
is the coefficient containing one- loop suppression factor and CP violating phases. The 
observed value of the baryon asymmetry, B ~ 10~ 10 , is then obtained without any fine 
tuning of parameters. 

Our findings have also important implications for the conjecture that ultra-high cos- 
mic rays, above the Greisen-Zatsepin-Kuzmin cut-off of the cosmic ray spectrum, may 
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be produced in decays of superheavy long-living particles 0, |10[. In order to produce 
cosmic rays of energies larger than about 10 13 GeV, the mass of the X-particles must 
be very large, Mx ~ 10 13 GeV and their lifetime tx cannot be much smaller than the 
age of the Universe, Tx > 10 10 yr. With the smallest value of the lifetime, the observed 
flux of ultra-high energy cosmic rays will be reproduced with a rather low density of 
X-particles, fix ~ 1CT 12 . It has been suggested that X-particles can be produced in the 
right amount by usual collisions and decay processes taking place during the reheating 



stage after inflation, if the reheat temperature never exceeded Mx |JT0[] . Again, assuming 
naively that that the maximum number density of a massive particle species X produced 
after inflation is suppressed by a factor of (M x /T R h) 3 ^ 2 exp(— M X /T RH ) with respect to 
the photon number density, one concludes that the reheat temperature Trh should be 
in the range 10 11 to 10 15 GeV 0. This is a rather high value and leads to the gravitino 
problem in generic supersymmetric models. This is one reason alternative production 



mechanisms of these superheavy X-particles have been proposed [11, 12, 13|. However, 
our analysis show that the situation is much more promising. Making use of Eq. (P2"D, 
the right amount of X-particles to explain the observed ultra-high energy cosmic rays is 
produced for 

( T «» ) ~ r^ 3/14 ( Mx \ (24) 

U0 10 GeW V200y U0 15 GeW' 1 ; 

where we have assumed (cr\v\) ~ M x 2 . Therefore, we conclude that particles as massive 
as 10 15 GeV may be generated during the reheating stage in abundances large enough 
to explain the ultra-high energy cosmic rays even if the reheat temperature satisfies the 
gravitino bound. 
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APPENDIX A: THERMALIZATION OF LIGHT DEGREES OF FREEDOM 

The form of the Boltzmann equations we use, e.g., Eq. @, assumes that the light 
particle decay products of the inflaton field are in local thermodynamic equilibrium 
(LTE). In this appendix we discuss this assumption, and the implications if it is not 
valid. 

Before discussing the validity of the assumption, it is useful to recall why the assump- 
tion was made. In the derivation of Eq. (fj), one starts with an equation for the rate 
of change of the X number density due to the process 77 — > XX with four-momentum 
conservation p 1 + p'^ = px + p'x- 



n x 



d% fd% r<fp x [<Pp' x 8 4 / , , n 

xfifa)M)\M\' xx + >- - • (Al) 



Here fi(p) is the phase-space density of particle species i with momentum pi, and 
^xx * s the sc i uare of the matrix element for the process 77 — > XX. With the 
assumption that the light particles are in LTE, the product of the light particle phase- 
space densities is / 7 (p 7 )/ 7 (p 7 ) = exp(— E^/T) exp(— E' /T). This last product is, of 
course, simply fx Q {Px)fx Q (p'x)i "which, after some rearrangement of Eq. (|A1|) , leads to 
a term for the creation of A's proportional to (^x^) ■ 

nx = {<r\v\) (n E x Q ) 2 + ■ ■ ■ . (A2) 

(for complete details, see Hj). 
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The factor (n x j in Eq. is present because not every light-particle collision 

has sufficient center-of-mass energy to create an X pair. If LTE is established with 
temperature T < Mx, the factor exp(— 2Mx/T) in (n x ®^J represents the fraction of the 
collisions with center-of-mass energy above threshold, i.e., with > 2Mx- 

A simple indication of whether thermalization occurs on a timescale shorter than 
the timescale for X production is the ratio of the cross section for the thermalization 
reactions to the cross section for X production. If the ratio is larger than unity, then 
thermalization of the light degrees of freedom is a good assumption. 

The process of X production involves a "hard" process, and the cross section will be 
acx/Mx, where ax was defined in Eq. (|]). In order to produce an equilibrium distribution 
from the original decay distribution it is necessary to change the number of particles. 
Therefore, the relevant cross section is the one for processes like 77 — > 777 .[] Although 
the thermalization reaction is higher order in perturbation theory, it is a "softer" process, 
and radiation of a soft photon has a large cross section. 

Without knowing the details of the interactions of the decay products, it is impossible 
to say with certainty how complete thermalization will be. But if the inflaton decay 
products have usual gauge interactions, the thermalization cross section will be larger 
than the X production cross section, and thermalization of the inflaton decay products 
is likely. 

Now let's explore the consequences if LTE of the light degrees of freedom is not 
obtained. If the light particles are not in LTE, then the factor n x lYi Eq. ( |A2|) could 
simply be replaced by the more general factor n^E > Mx)-f\ Now let's make the extreme 

assumption that the light degrees of freedom never interact before X production, and 

4 Recall that 7 represents a light particle, not just a photon, so 7 may carry electric charge, color 
charge, etc. 

5 Of course in this case would not be a thermal average, but an average over the actual phase- 

space density. 
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that they have the original (redshifted) momentum with which they were created in 
inflaton decay Assume this original momentum is M^/rj. If M^/r] is greater than 
Mx, then the X production rate actually will be larger than the equilibrium rate since 
n^{E > Mx) > n x , while if M^/rj is less than Mx, then the X production rate will be 
zero! The most reasonable assumption is that even if there is a large multiplicity in 
decay, a fair fraction of M^ is carried by a few leading particles. So the effective value 
of 77 is probably not too large. 

The above analysis leads us to the conclusion that thermalization of the light degrees 
of freedom is likely unless the inflaton decay products themselves are very weakly coupled 
to everything. Even if thermalization does not occur, production of massive particles 
during reheating is not much different than our simple model suggests. 
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